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Abstract. - The Harper-Hofstadter problem for two-dimensional electron gas with Rashba 
spin-orbit coupling subject to periodic potential and perpendicular magnetic field is studied 
analytically and numerically. The butterfly-like energy spectrum, spinor wave functions as 
well as the spin density and average spin polarization are calculated for actual parameters of 
semiconductor structure. Our calculations show that in two-dimensional electron gas subject 
to periodic potential and uniform magnetic field the effects of energy spectrum splitting caused 
by large spin-orbit Rashba coupling can be observed experimentally. 



Introduction. - The problem of quantum states of an electron subject to both periodic 
potential and homogeneous magnetic field remains to be actual for several last decades. There 
are various approaches and models for periodic potential and different approximations which 
are used for investigations of these states (see, ref. [1]). However, the spin-orbit interaction is 
usually excluded from the models as well as the Zeeman term. Such approach can be justified 
as long as the amplitude of the periodic potential is big enough to make the Landau level 
splitting much greater than the typical spin-orbit coupling energy and the Zeeman term. At 
the same time, under the realistic experimental conditions with 2D electron gas subject to the 
potential of lateral superlattice the potential amplitude Vb can be of the same order as the 
spin-orbit (SO) Rashba coupling. For example, in recent papers [2,3] where a pioneering step 
into experimental observation of 2D magnetic Bloch states has been done, the magnitude of 
Vo was around 1-5 meV. In semiconductors structures with large SO coupling [4] the typical 
splitting energy can be of the same order as Vo. 

In this paper we study magnetic Bloch states of 2D electrons subject to both periodic 
potential of a lateral superlattice and perpendicular magnetic field under the conditions when 
the SO coupling and Zeeman term should be taken into consideration. The energy band 
structure is calculated in the magnetic Brillouin zone (MBZ) and the magnetic Bloch states 
are constructed. Also, the spin density distribution in the elementary cell is obtained and the 
average spin polarization in a state with given quasimomentum are calculated. It is shown that 
the consideration of spin-orbit coupling is necessary for interpretation of realistic experiments. 
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Hamiltonian matrix structure, energy band spectrum and wave junctions. - We consider 
2D electron gas with spin-orbit Rashba coupling in a potential V(x, y) which is periodic in 
plane with the period a, and in a uniform magnetic field H perpendicular to the plane of the 
electrons. The correspondent one-electron Hamiltonian has the following form: 

H = H + V(x,y), (1) 

where V(x, y) = V(x + na, y + ma) is a periodic potential, 



H = (p - eA/c) /2m* + - \a x (p y - eA y /c) - a y p x j - g[i B Ha z 

is the Rashba Hamiltonian of an electron in uniform magnetic field [5,6]. Here, p x . y are the 
momentum operator components, m* is the electron effective mass, a are the Pauli matrices, a 
is the parameter of the SO coupling, g is the Zeeman factor, and \xb is the Bohr magneton. We 
use the Landau gauge in which the vector potential has the form A = (0, Hx, 0) and consider 
the potential V(x,y) = Vo(cos(2Trx/a) + cos(2ny / a)) . The quantum states structure of the 
system under consideration depends crucially on the parameter $/$o = p/q = \e\Ha 2 /2ttUc 
(p and q are prime integers) which is the number of flux quanta per unit cell, and <J> is the 
flux quanta. 

One can express the eigenfunction of Hamiltonian (1) as a set of Landau wave functions 
in the presence of the SO coupling [7]. If we take p/q to be the rational number, the two- 
component magnetic Bloch function, analogous to the one-component magnetic Bloch function 
of the paper [1] , can be written in the form 

^k(x,y) — ( lk '- a ''^-' J — e ik y y e ik x (iqa+nqa,/p) e 2niy(ep+n)/a x 
' " ' S=ln=li=-oo 



^2k(x,y) 



(2) 



where the spinors ^ k = exp(^) (j^), tfU = ^gf {^[S^) 
■4,- = ex P (ifc ^ } ( , ] correspond to "+" and "-" branches of the spectrum of 



Snfl[ V^s \-DsMtU*] 

the Hamiltonian H [7], D s = {V2Sa/l H )/(E^ + \j {E^) 2 + 2Sa 2 /l 2 H ), <j> s [€\ is the simple 

harmonic oscillator functions, Ih = ch/\e\H is the magnetic length, Eq = hui c /2 + g/J-sEt 
and uj c = \e\H/m*c is the cyclotron frequency. Here quantum numbers S = 1, 2, 3, . . . charac- 
terize the pair of "+" and "-" states Ef = SHlj c ± ^ (E+) 2 + 2Sa 2 /l 2 H of the unperturbed 

Hamiltonian Ho, & n = (x — xo — iqa — nqa/p)/lh, xo — chk y /\e\H. 

Note that the spinor wave function (2) is the eigenfunction of both the Hamiltonian (1) 
and operator of magnetic translation and therefore it have to obey the following Bloch-Peierls 
conditions 

*k(a; + qa, y + a) = ^ k (x, y) enp(ik x qa) exp(ik y a) exp(27ny/a), 
where k is quasimomentum defined in the MBZ 



—n/qa < k x < n/qa, —n/a < k y < n/a. 
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So, the magnetic Brillouin zone is the same as for charged spinless particle. At the limit of 
high magnetic fields the functions V'snflc an d ^Snik. are proportional to eigenspinors of Pauli 

operator a z : {]^J an d ^ ^ 

Substituting eq.(2) in the Schrodinger equation H^f = we come to the infinite system 
of linear equations for coefficients ^4s„(k) and i?5„(k). In the case when periodic potential 
amplitude and spin-orbit coupling energy have the same order and the inequality AEso — 
Vq < hu c takes place the system of linear equations can be reduced to the system of infinite 
number of uncoupled groups of 2p equations. Each group describes the magnetic Bloch states 
formed from the states with energies Eg and E$ +1 of the single Landau level split by SO 
interaction. 

In this approximation the system of 2p linear equations attached to the quantum number 
S is defined by the following Hamiltonian matrix 



/Gi 

M* 



rrSSr _ 
n nnl ~ 



M 
G 2 
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M* 
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M* G p J* ... 
... J C/i N ... 
... N* U 2 N 



± p 

N* 
... 











N* 



(3) 



N* U V J 



p " " • • • " *■ 
which has the block structure. In eq.(3) the matrix elements are defined as follows 

G n = E+ + V c-^/ 2 Pcos(27Tx /a + 27rnq/p)[D 2 s L° s _ 1 (wq/p) + L° (Trg/p)]— Lj 
M =Yl eMik x qa/p)e-^' 2 nD 2 S L%_ 1 {-Kq/p) + L§ (7rg/p)] _J_ 

U n = E- +1 + V a c-^/ 2 Pcos(27rx /a + 2nnq/p)[D 2 s+1 L° s+1 (nq/p) + L° s (nq/p)} 1 

1 + u s+ 

N = ^e ik ^e-^[D 2 s+1 L s+1 (nq/p) + L^q/p)] ^ 1 ^ , 



(4) 



F n = V e~ v * sin(27r^ + 2tt^) 
a p 



(1 + D|)(1 + D§ +1 )' 



V2 y 2 P 

V2 



(1 + D|)(1 + £>| +1 )' 
1 



(i + D 2 )(i + D 2 s+1 y 



Here if (z) is the Laguerre polynomial. 

This system has a form of generalized Harper's equation. Matrix elements G n and M in 
the left upper block of eq.(3) characterize the interaction between different states of the "+" 
branch and elements U n and N in right low block are due to the interaction between states 
of the "— " branch. The right upper and left low blocks with elements F n , J and T describe 
the transition between states of the neighbor "+" and "— " branches with index S. Matrix 
elements in these four blocks are periodic in n with period p. So, each block has the size pxp. 
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Fig. 1 - Position of energy subbands versus magnetic flux numbers at the parameters indicated in 
the text. The big dots show "+" and "— " Landau level positions in the absence of periodic potential. 
The insert shows 2p energy subbands for p/q — 3/1, 4/1, 5/1. 



Numerical calculations of energy spectrum and wave functions were carried out for the 
actual parameters which approximately correspond to the structure used in the experiment. 
In our calculations we used the following parameters of 2D electron gas: Vo = 1 meV, effective 
mass to* = 0.05 mo, g — ginAs = —2.0, a = 5 • 10 -11 eV ■ to. The lattice period was chosen as 
a = 60 Tito and Vb = 1 meV. Under these parameters the corresponding spectrum structure 
in magnetic field becomes well resolved and SO splitting is larger than Zeeman splitting at 
p/q < 5. The chosen parameters approximately correspond to the InAs quantum well where 
the Rashba SO constant can reach the maximal value 5 • 10 -11 eV ■ to due to electron wave 
function penetration into barrier layer [4]. 

In fig.l we show the electron energy levels for the state k = versus the magnetic flux 
number. The big dots define the position of Landau levels Eg in the absence of periodic 
potential. The insert demonstrates the magnetic subbands at p/q equal to 3/1, 4/1 and 5/1 
for k located in the magnetic Brillouin zone. The upper arrow marks the lowest magnetic 
subband for which the dispersion law Ei(k) is plotted in fig. 2. Note that the spectrum of 
magnetic subbands has the symmetry of the C^v group, as it should be. As one can see from 
fig.l the periodic potential forms the spectra which resemble the Hofstadter butterflies at the 
region 1 < p/q < 2. For p/q > 3 the magnetic subbands split by SO coupling and periodic 
potential are not overlapped. So, at magnetic fields corresponding to p/q > 3 the unperturbed 
levels (fat dots in fig.l) are grouped into pairs and, therefore, our two-branch approximation 
is appropriate. 

To compare the effects of level splitting due to the periodic potential as well as Zeeman 
effect and SO coupling we calculated the energy spectra formed from the lowest pair of Landau 
levels (see fig. 3). In fig. 3a we show the Hofstadter-like spectrum in the absence of Zeeman and 
SO interactions. Here, all levels are two-fold degenerate. At the presence of Zeeman interaction 
(fig. 3b) the degeneracy is lifted and the correspondent spectrum is the superposition of two 
Hofstadter-like butterflies. Fig. 3c demonstrates the splitting due to the SO Rashba coupling. 
So, under the parameters and magnetic fields indicated above the SO interaction is large 
enough to form two non-overlapped and non-degenerate groups of magnetic subbands. 

We have calculated also the electron density |\E'k(a ; j y)\ 2 — \^ , y)\ 2 + \^2^{ x tU)\ 2 at 
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Fig. 2 - Electron energy versus quasimomentum in magnetic Brillouin zone. 



k = for the fourth magnetic subband labelled by low arrow at the insert in fig.l. At the 
center of MBZ it has the symmetry of the C^v group. Note, that for the state k ^ the 
electron density does not possess the symmetry of the lattice. 

Spin density and average spin polarization. In the 2D electron gas with spin-orbit 
coupling at the presence of the periodic potential and perpendicular magnetic field the non- 
trivial spin structure appears in each of magnetic energy bands. Such a structure at a certain 
k can be characterized by the vector field of the spin density in coordinate space 

S ik (x,y) = 4-+(r)o;* k (r), (i = x,y). 

In addition, spin states can be characterized by the average spin distribution in k space 

5i(k) =< S ik (x,y) >=< #£(r)oi# k (r) >, 

where the brackets mean the integration on x and y over the magnetic unit cell. 

In fig. 5 we show the spin density distribution calculated in the r-space for the square 
superlattice with parameters indicated above. All results correspond to the fourth magnetic 
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Fig. 3 - Energy spectra in the absence of SO Rashba coupling and Zeeman effect (a), in the presence 
of Zeeman splitting only (b), and when all of the contributions to level splitting are present (c). The 
parameters are the same as in fig.l . 
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Fig. 4 - The electron density \^ k (x,y)\ 2 = \^i k (x,y)\ 2 + \^ 2 k(x,y)\ 2 at k = for the fourth magnetic 
subband. 



energy band shown at the insert in fig.l for the state k = 0. The S x k=o(%, y) and S y k = o(x, y) 
densities have the C s symmetry, but the S z k = o(x,y) component has the C^v symmetry. The 
integration of the spin densities over the magnetic unit cell leads to the zero values of S x (k) 
and S y (k) at k = (see fig. 6a). The integration of S z -k = o(x, y) over the magnetic unit cell 
results to the positive value of the average spin projection along magnetic field. 

In fig. 6 we show the average spin polarization in the MBZ for the fourth band at the 
magnetic flux p/q = 3/1 (sec fig.l). Fig. 6a demonstrates the distribution of the components 
5a; (k) and S y (k) and fig. 6b shows the iS'z(k) component. The distribution in fig. 6a has the 
vortex structure with the C2 symmetry. The vortex centers are located at the center and at 
the corners of MBZ. These two vortices have opposite directions of rotation. In these points 
the average lateral spin components are smaller than perpendicular component S^(k). For 
the state k = such a behavior is due to the small value of SO coupling in comparison 
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Fig. 5 - Electron spin densities in a unit cell of the lateral superlattice. 



V. Ya. Demikhovskii and A. A. Perov: Harper-Hofstadter problem for 2D electron gas with Ic-linear Rashba spin-( 



7I/qa 


< i i i i i i ■ 

- < « i t / /■ / * 

-^^///^^-- 

' ' / / / / / / 

-'■/////y- 

' ■* , s / / s ^ 

'-'// / ./ / / -~ ~ 
-'////// ^ - 
' ' i t 1 I J I ' " 


\ \ \ \ \ \ \ \ • 
\ \ \ \ \ \ . 
xN\\\\\n, _ 
~-XX\X\\v-> _ 

-,xX\X\\-x-._ 
xXX\\\\x-,_ 

xX\\\\\N-, _ 

x\\\\\\\n. 
X \ \ \ \ \ \ - 

.III 


-,t/qa S z (k).a.u. 

/a 

-It/a — - — o 

— ^x^n/qa 




1 i l 1 1 1 l 

s * * \ \ \ \ \ v 

-->^\\\\\N~- 
-~-x v %\XXX~ 

- i X V \ v \ •< ■ 


_///////,< 

X /////// , 

-X/// / / J 

_- x X / / / / / / 

.. 

__^XXX///^^ 
_ x ■/ / / / / / x, 
-,///////, 
,////([/, 

t * ; i i < ! 


ft/qa ' .-. ■ ■ ' . ' ' . . ■ 



(a) <t» 

Fig. 6 - (a) Distribution of average spin polarization S x (k) and S y (k) in the fourth magnetic subband 
attached to the two lowest neighbor branches, (b) 5* z (k) in the fourth magnetic energy subband. 



with Zeeman effect. At the corners of MBZ the SO coupling is effectively smaller than the 
spin splitting in a magnetic field as a result of state mixing at the magnetic Brillouin zone 
edges. Here, in accordance with normalizing condition the |.!5 z (k)| component has the local 
maxima. We would like to note here that these spin vortices define the Berry phases and, 
respectively, the spin Hall conductance. Due to the same fact of state mixing by periodic 
potential the SO coupling effectively vanishes at the points where the edges of MBZ cross the 
k x — and k y = lines. Finally, note that the total magnetic moment of 2D electron gas 
which completely occupies this zone is positive. 

Conclusion. - We have considered the effects of spin-orbit coupling on the quantum states 
of two-dimensional electron gas at the presence of periodic potential and uniform magnetic 
field. The Hofstadter butterfly-like energy band spectrum and electron density in the unit cell 
was analyzed. We demonstrate that the Zeeman interaction and the spin-orbit coupling form 
a complex and intricate spin periodic structure in the field of the lateral supcrlattice. The 
electron spin densities S x k(x,y), Sy^Xjy) and S z ^{x,y) as well as average spin polarizations 
Sa;(k), ^(k) and S z {k) for the states with different quasimomentum k defined in magnetic 
Brillouin zone were calculated and the role of three contribution to the spin orientation was 
investigated. We think that our results will be useful for analyzing the data of real experiments 
on the hunting of Hofstadter butterfly in two-dimensional electron gas. 
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